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nical Monitor  was  Mr.  George  Matsumura,  DAEN-MCE-A. 
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DYNAMIC  RESPONSE  OF  REINFORCED 
CONCRETE  STRUCTURES 


1 INTRODUCTION 


Background 

Reinforced  concrete  (RC)  construction  is  one  of  the  major  types 
of  construction  used  throughout  the  world.  At  present,  the  design  and 
analysis  of  most  reinforced  concrete  structures  are  based  on  linear 
elastic  procedures,  even  though  some  of  these  structures  will  be 
subjected  to  loads  in  the  inelastic  range.  Structures  in  severe  earth- 
quake zones,  for  example,  may  undergo  many  cycles  of  inelastic  deforma- 
tions during  their  service  lives.  Thus,  there  is  a need  for  refined 
cyclic  inelastic  analyses  of  structures  in  order  to  develop  a more 
economical,  realistic  approach  and  safer  design  procedures.  Such 
analyses  require  knowledge  of  the  behavior  of  structural  materials 
under  arbitrary  complex  cyclic  loading. 


Obj ective 


The  objective  of  this  study  is  to  develop  criteria  for  determining 
damping  coefficients  and  ductility  factors  in  RC  structures  based  on 
their  energy  dissipation  characteristics. 

The  study  is  divided  into  three  phases: 

1.  Development  of  a uniaxial  model  for  the  cyclic  inelastic 
response  of  reinforced  concrete. 

2.  Development  of  a beam-column  finite  element  and  a computer  pro- 
gram for  determining  inelastic  response  for  RC  structures  under  seismic 
loads . 

3.  Extension  of  the  computer  program  for  estimating  energy  dissipa- 
tion and  equivalent  damping  coefficients  and  ductility  factors. 

The  first  phase  of  study  has  been  completed  and  is  described  in 
CERL  Technical  Report  M-180.'  This  report  covers  the  work  accomplished 
in  the  second  phase. 


'Sharma,  S.  K.  and  Bhattacharyya , R.  K.  , An  Analytical  Model  for  Uniaxial 
Cyolia  Inelastic  Behavior  of  Reinforced  Concrete ^ Technical  Report 
M-180/ADA024910  (U.S.  Army  Construction  Engineering  Research  Laboratory 
[CERL],  1976). 


5 


Beam-column  and  plane  stress  elements  are  developed  using  cubic 
Hermite  polynomials.  A constant  acceleration  method  (a  special  case 
of  Newmark's  6 method)  is  used  in  the  solution  procedure.  The  method 
assures  a constant  acceleration  between  time  steps.  The  elements  are 
incorporated  in  a commercially  available  ■nnuter  program,  DRAIN-2D, 
and  a two-story,  one-bay  reinforced  concrtce  structure  is  analyzed  as 
an  example  problem. 


Scope 


The  RC  structures  considered  in  this  report  are  plane  structures 
(two-dimensional)  with  inelastic  beams  and  columns,  and  elastic  shear 
walls.  The  shear  walls  may,  however,  fail  in  a brittle  fracture  mode. 
Structural  elements  which  fail  primarily  by  brittle  fracture  show  a 
relatively  linear  elastic  behavior  up  to  the  failure  load.  Hence 
there  would  be  no  significant  increase  in  accuracy,  but  possibly  a 
noticeable  increase  in  computational  time  and  cost,  if  the  inelasticity 
of  shear  walls  was  also  considered.  The  structures  are  idealized  by 
beam-column  (for  beams  and  columns)  and  plane  stress  (for  shear  walls) 
finite  elements. 


Mode  of  Technology  Transfer 

The  criteria  developed  in  this  study  are  to  be  incorporated  into 
a seismic  design  manual,  TM  5-809-10,  Seismic  Designs  for  Buildings. 


2 BEAM-COLUKH  AND  PLANE  STRESS  FINITE  ELEMENTS 


The  existing  finite  element  computer  programs  for  RC  structures 
(e.g.,  SAKE,  STRUDL,  PRIEST,  CONCRETE,  etc.)  use  elastic  beam-column 
elements  with  nonlinear  moment-rotation  springs  at  the  ends.  The 
properties  of  these  nonlinear  springs  are  chosen  to  simulate  stiffness 
changes  caused  by  cracking  of  the  concrete,  yielding  of  the  rein- 
forcement, and  stress  reversals.^  Implicit  is  the  assumption  that 
inelastic  deformations  are  concentrated  only  at  the  ends  of  beams  and 
columns.  However,  this  is  too  restrictive.  In  a real  RC  structure 
inelastic  deformations  can  occur  throughout  the  member.  Thus,  energy 
dissipation  estimates  based  on  these  beam-column  elements  are  not 
expected  to  be  reliable.  Consequently,  a new  beam-column  finite 
element  was  developed  which  permits  cyclic  inelastic  deformations  at 
all  points  in  the  element. 

The  following  assumptions  are  employed  for  the  structure,  pre- 
scribed forces,  reinforcement,  and  beam-column  and  plane  stress  finite 
elements. 

1.  Only  two-dimensional  RC  structures  are  considered. 

2.  Ground  motion  is  parallel  to  the  plane  of  the  structure. 

3.  The  structures  are  fixed  at  their  bases  on  infinitely  rigid 
foundations. 

4.  Mass  of  the  floors  is  assumed  to  be  concentrated  at  the 
member  joints  (nodal  points). 

5.  Displacements  and  strains  are  assumed  to  be  small. 

6.  Shear  deformations  in  beams  and  columns  are  ignored. 

7.  There  is  no  slip  between  the  concrete  and  the  reinforcing 
steel . 

8.  Kirchhoff's  assumption,  which  states  that  normals  to  the 
centroidal  axis  remain  straight  and  normal,  is  used. 

9.  Reinforcing  steel  is  "smeared"  in  the  concrete  to  obtain  an 
equivalent  tangent  modulus. 


^Otani , S. , SAKE,  A Computer  Program  for  Inelastic  Response  of  R/C 
Frames  to  Earthquakes,  Civil  Engineers  Studies,  Structural  Research 
Series  No.  413  (University  of  Illinois  at  Urbana-Champaign , November 
1974). 
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Beam-Column  Element 

Figure  1 shows  a beam-column  finite  element.  In  the  local 
coordinate  system,  the  centroidal  axis  of  the  beam  is  on  the  x-axis 
from  X = 0 to  X = L.  Displacements  along  x and  y axes  at  any  point 
X on  the  centroidal  axis  are  denoted  by  u(x)  and  v(x),  respectively. 
These  displacement  components  are  represented  as 


u = U.|  + 

V 


[Eq  1] 


ViTii  + V-n2  + V2n3  + V'n 


where 


Up  V-j  = displacement  components  at  node  1 

U^,  ^2  ~ displacement  components  at  node  2 

VI,  Vp  = derivatives  with  respect  to  x of  v at  nodes  1 and  2 
respectively 

tjj  = a linear  function  of  x defined  in  Appendix  A 

rii  ,ri2.n3»n4  = cubic  Hermite  polynomials  of  x defined  in  Appendix  A. 

Axial  strain  at  any  point  (x,  y,  z)  of  the  beam  is  given  by 
Sloane. ^ 


e^(x,y,z)  = E^(x,y,0)  - u'  - yv" 


[Eq  2] 


where 

= strain  at  any  point  in  the  x direction 
V',V"  = derivatives  with  respect  to  x 
therefore 


6e^  = 6u'  - y6v" 

where  (S  denotes  a virtual  increment. 


[Eq  3] 


^Sloane,  A.,  Mechanias  of  Materials  (Dover  Publications,  Inc.,  1952). 
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Therefore,  following  the  notation  of  Zienkiewicz'*  yields 

[B]  = [0  - yn'i  - yn^  <1^'  - yn;;  - yn;;]  [Eq  « 

The  tangential  stiffness  matrix  for  the  beam-column  element  ma. 
now  be  written  as^ 


[k^]  = [B^  [E^][B]dv  [Eq  ’ 

‘‘Zienkiewicz,  0.  C.,  The  Finite  Element  Method  in  Engineering  ■ 

(McGraw-Hill,  1971). 

^Zienkiewicz . 


where 


[k-j.]  = tangential  stiffness  matrix 
[B]^  = transpose  of  [B]  defined  by  Eq  6 
= tangential  modulus 
V = volume  of  the  beam-column  element. 


Assuming  the  breadth  of  the  beam  to  be  b at  a distance  y from  the 
centroidal  axis,  Eq  7 may  be  written  as  ^ 


[k^]  = 


0 Yi 


[B]  [E.p][B]  bydydx 


[Eq  8] 


The  integral  of  Eq  8 is  evaluated  by  numerical  integration  using 
Simpson's  rule.  A 5 x 5 matrix  of  integration  points  is  chosen  across 
the  depth  and  length  of  the  beam. 

The  tangent  modulus  at  an  integration  point  (x^,y^)  is  given  by 


Et  = u(Xi.yi)  E^ 

{1  - u(x.,y.)}  E^ 


[Eq  9] 


where 

y(x-,y-)  = reinforcement  ratio  (ratio  of  cross  sections  of  rein- 
''  forcement  area  to  the  total  area)  at  (x^.  ,y^. ) 

E^  = tangent  modulus  of  the  reinforcing  steel  at  (x^  ,y^ ) 

E^  = tangent  modulus  of  the  concrete  at  (x.,y.). 

The  tangent  moduli  of  the  concrete  and  the  reinforcing  steel  are 
obtained  by  their  cyclic  inelastic  models  as  described  in  Sharma  and 
Bhattacharyya.® 


'’Sharma,  S.  K.  and  Bhattacharyya , R.  K.,  An  Analytiail  Model  for  Uni^ixial 
Cyolia  Inelastic  Behavior  of  Reinforced  Concrete,  Technical  Report 
M-180/ADA024910  (CERL,  1976). 
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The  tangential  stiffness  matrix  in  global  coordinate  system 
{X,Y,Z)  is  given  as^ 


[K^]  = [r]^[4][r] 


[Eq  10] 


where 

[Kj]  = tangential  stiffness  matrix  in  global  coordinate  system 

[r]  = transformation  matrix  from  local  to  global  coordinate 
system 

[r]^  = transpose  of  [r]. 


The  transformation  matrix  [r]  is  given  in  Appendix  B. 


Plane  Stress  Element 


A plane  stress  finite  element  is  used  to  discretize  shear  walls. 
Several  highly  accurate  plane  stress  elements  are  available  in  the 
literature  and  one  of  these  elements®  has  been  selected.  This  plane 
stress  element  (Figure  2)  is  a triangular  element  with  three  nodes. 
Displacement  components  u and  v within  the  element  are  written  as 


u = a-j  + a2X  + a^y 

V = 04  + agX  + a^y 


[Eq  11] 


where  (x,y)  are  local  coordinates  and  •••>  “5  generalized 

nodal  variables.  Substituting  nodal  coordinates  (X.,y.),  i = 1,  2,  3 
in  Eq  11,  a,,  02,  ...,  ag  can  be  expressed  in  terms^of  u^v,  ,U2,V2,u- 
and  v^,  where  (u.,v.)  denote  the  displacement  components  at  node  t, 
i = 1,  2,  3.  This  results  in  the  following  strain-displacement  relation 


^Zienkiewicz , 0.  C.,  The  Finite  Element  Method  in  Engineering  Soienae 

(McGraw-Hill,  1971). 

®Zienkiewicz . 

®Zienkiewicz. 
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and 


A = area  of  the  triangle  with  nodes  1,  2,  and  3. 
Consequently,  the  strain-displacement  matrix  [B]  is  given 


0 

c. 


The  tangent  stiffness  matrix  is  now  written  as^® 


[K^] 


[B]^ [D^][B]  dx  dy 
A 


where  [B]  is  given  by  Eq  13,  [B]^  is  the  transpose  of  [B],  and 
is  the  plane  stress  elasticity  matrix; 


[D^] 


V 0 
1 0 

0 (l-v)/2 


‘ ®Zienkiewicz , 0.  C.,  The  Finite  Element  Method  in  Engineering 

(McGraw-Hill,  1971). 


[Eq  12] 


by 

[Eq  13] 

[Eq  14] 

[D^] 

Science 
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where 

Ej  = Young’s  modulus  of  elasticity 
V = Poisson's  ratio. 


The  plane  stress  element  is  assumed  to  be  elastic,  but  may  fail  in 
tension  (brittle  failure).  The  failure  of  an  element  in  tension  is 
modeled  by  setting  its  elastic  modulus  E-r  to  zero  once  the  principal 
tensile  stress  in  the  element  reaches  a prescribed  value. 


3 SOLUTION  PROCEDURE 


At  any  instant  of  time  t,  the  incremental  equation  of  equilibrium 
may  be  written  as** 

[M]  {dr}  + [K^]  {dr}=.  {dP}  + {C^}  [Eq  16] 

where 

[M]  = lumped  mass  matrix  of  the  structure 

{df}  = increment  of  nodal  acceleration  matrix 

[K^]  = structural  tangential  stiffness  matrix  at  time  t 

{dP}  = increment  of  externally  applied  force  matrix 

{C  } = load  correction  matrix 
r 

{dr}  = increment  of  nodal  displacement  matrix 


Eq  16  is  solved  by  using  a constant  acceleration  method  (one  of  the 
special  cases  of  the  Newmark  $ method).  This  method  assumes  a constant 
acceleration  between  the  time  steps  t and  t+dt. 

Therefore, 


{f}  = 1/2  {(r^ 


* 


[Eq  17] 


where 


{f}  = assumed  constant  acceleration  between  t and  t+dt 
{f^}  = acceleration  at  time  t 
{f’t^dt^  = acceleration  at  time  t+dt 


“McNamara,  J,  F.  and  Sharma,  S.  K.,  An  Analytical  Model  for  Determining 
Energy  Dissipation  in  Dynamically  Loaded  Structures,  Technical  Report 
M-165/ADA017040  (CERL,  1975). 
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Eq  17  is  rewritten  as 


where 


= h * f } 


= ^Vdt  - ■•‘t^ 


[Eq  18] 


[Eq  19] 


Integration  of  Eq  18  yields 

{dr)  • |-2rj  t 


{dr} 


4r 


4dr 


(dt)' 


[Eq  20] 


where 


{f^}  = nodal  velocity  vector  at  time  t 


Substituting  Eq  20  into  Eq  16  yields 


(dt) 


2 [M]  + [K^] 


{dr}  = {dP}  + {C^}  + [M]  |2f^  + [Eq  21] 


Eq  21  is  solved  for  {dr},  which  upon  substitution  in  Eq  20  yields 
{dr}  and  {df}.  Thus,  solving  Eqs  21  and  20  successively  for  each  time 
step  yields  time-history  solutions  for  displacement,  velocity,  and 
acceleration  of  the  structure. 
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4 RESULTS 


DRAIN-2D 


The  beam-column  and  plane  stress  finite  elements  described  in  this 
report  have  been  incorporated  in  a commercially  available  computer 
program,  DRAIN-2D.*^  Other  finite  elements  in  the  element  library 
DRAIN-2D  are  (1)  an  elasto-plastic  truss-element,  (2)  a beam  element 
with  plastic  hinges  at  the  ends,  (3)  a semi-rigid  connection  element, 
and  (4)  a shear  panel  element  which  has  shear  stiffness  only.  This  pro- 
gram can  be  used  for  both  static  and  dynamic  (linear  or  nonlinear) 
analyses  of  plane  frame  structures.  However,  dynamic  analyses  are  re- 
stricted to  earthquake  excitation  in  which  all  support  points  are 
assumed  to  move  in  phase.  The  earthquake  excitation  is  defined  by  time 
histories  of  ground  acceleration  which  may  have  both  horizontal  and 
vertical  components.  A two-story,  one-bay  reinforced  concrete  structure 
(Figure  3a)  is  analyzed  by  DRAIN-2D  in  the  following  example. 


Example  Problem 

Figure  3b  shows  a typical  cross-section  of  beams  and  columns  of 
the  structure.  Hatched  areas  of  the  cross-section  are  steel  reinforced. 
Table  1 gives  the  dimensions  of  these  beams  and  columns  and  reinforce- 
ment ratios  of  the  hatched  areas.  Dead  weights  on  each  story  of  the 
structure  are  50  kips  (22,680  kg). 

The  structure  is  subjected  to  a north-south  component  of  the  El 
Centro  earthquake.  The  acceleration-time  history  of  the  earthquake 
(Figure  4)  is  used  as  the  base  acceleration  for  the  structure. 

Stress-strain  response  of  the  concrete  and  the  reinforcing  steel 
is  predicted  by  their  uniaxial  cyclic  models  described  in  Sharma  and 
Bhattacharyya.  These  models  have  been  programmed  into  DRAIN-2D  to 
obtain  tangent  moduli  Eg  and  E^  (see  Eq  9)  at  each  time  step  of  the 
numerical  integration  algorithm.  A time  step  of  0.01  seconds  is  used 
for  numerical  integration,  and  results  are  printed  out  after  every 
10  steps. 

The  displacement-time  response  of  node  5 (Figure  3a)  is  plotted  in 
Figure  5.  The  maximum  and  minimum  displacements  for  this  node,  which 
occur  at  2.2  sec.  and  2.7  sec.,  respectively,  are  1.14  in.  (29  imi)  and 
0.67  in.  (17  mm).  The  maximum  bending  moment  is  1414.48  kips-in. 

(16  297  000  kg-mm)  which  occurs  at  node  4 of  member  5 (Figure  3a)  at 
time  2.2  sec. 

^^Kanaan,  A.  E.  and  G.  H.  Powell,  DRAIN-2D,  A General  Purpose  Corrqmter 
Program  for  Dynamic  Analysis  Inelastic  Plane  Structures  (Earthquake 
Engineering  Research  Center,  College  of  Engineering,  University  of 
California,  Berkeley,  April  1973). 
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Figure  3.  Example  problem. 


Table  1 


Properties  of  the  Frame  Members 


Member 

Length,  L 
in.  (mm) 

Breadth,  b 
in.  (mm) 

Depth,  d 
in.  (mm) 

Reinforcement 
Ratio,  p 

1 

240.0 

12.0 

24.0 

0.04 

(609.6) 

(305) 

(610) 

2 

240.0 

10.0 

20.0 

0.04 

(609.6) 

(254) 

(508) 

3 

120-0 

12.0 

24.0 

0.04 

(304.8) 

(305) 

(610) 

4 

120.0 

10.0 

20.0 

0.04 

(304.8) 

(254) 

(508) 

5 

120.0 

12.0 

24.0 

0.04 

(304.8) 

(305) 

(610) 

6 

120.0 

10.0 

20.0 

0.04 

(304.8) 

(254) 

(508) 
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5 CONCLUSIONS 


Beam-column  and  plane  stress  finite  elements  have  been  described 
for  an  inelastic  analysis  of  plane  RC  structures  under  earthquake- 
type  ground  motion.  Material  nonlinearities  in  the  beam-column  finite 
element  are  taken  into  account  by  considering  cyclic  inelastic  de-  i 

formations  throughout  the  element.  The  plane  stress  finite  element 
allows  for  cracking  of  the  element  in  tension.  These  elements  have 
been  incorporated  in  the  DRAIN-2D  computer  program,  which  is  both 
flexible  and  efficient.  Preliminary  results  for  simple  structures 
show  that,  with  the  addition  of  new  finite  elements  described  in  this 
report,  this  program  will  be  very  useful  for  practical  investigations 
of  RC  structures. 

In  the  third  phase  of  this  research  project,  this  computer  pro- 
gram is  to  be  extended  by  adding  new  subroutines  for  the  calculation 
of  dissipated  energy  from  increments  of  strains  and  stresses  in  the  ; 

elements.  Numerical  estimates  of  dissipated  energy  will  then  be  used 
in  determining  damping  coefficients  and  ductility  factors  for  RC 
structures. 

The  capability  of  RC  structures  to  absorb  energy  through  material 
damping  and  ductility  contributes  significantly  to  preventing 
catastrophic  failure.  This  work  will  ultimately  provide  a scientific 
basis  for  estimating  damping  and/or  ductility  properties  for  the 

structural  elements  considered.  Other  elements  may  also  be  considered  i 

in  the  future,  using  similar  analysis  methods.  Accounting  for  energy 

absorbtion  with  sufficient  accuracy  in  this  manner  is  recommended  and 

will  permit  cost  savings  by  minimizing  design  stress  levels  or  load 

factors. 
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APPENDIX  A: 


FUNCTIONS  \p,  n-] » n2»  n^,  and 

The  functions  tp,  n-]  ♦ n2>  n3,  and  of  Eq  1 are  defined  as  follows 

4;  = (U^  - U^)  x/L 

q^  = + C2X  + 

^2  “ Cg  + CgX  + c^x  + CgX 

'^3  " *^9  ^11^  ^^12^^ 

■^4  ^ '^IS  ^15^  ^ ^16^^ 

where 


C-j  1 Cg  0 


C2  = 

0 

^10  " 

0 

"3  = 

-3/L^ 

3/L^ 

"4  = 

2/L^ 

C^2  = 

-2/L 

'^s  " 

0 

^13  ^ 

0 

1 

^14  " 

0 

-2/L 

^15  " 

-1/L 

^8  " 

1/L^ 

^16  = 

1/L^ 

1 
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APPENDIX  B; 


TRANSFORMATION  MATRIX  [r] 

The  transformation  matrix  [r]  is  given  by 

cos  6 sin  000  00 

-sin  0 cos  000  00 

0 0 1 0 0 0 

[r]  = 

0 00  cos  0 sin  0 0 

0 00  -sin  0 cos  0 0 

0 0 0 0 0 1 

where  0 is  defined  in  Figure  1. 
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